oo 

o 

O 
(N 

D 



in: 



> 
00 

o 

00 

o 



X 



Vector Meson Dominance as the first order 
of a sequence of Pade Approximants 

J.J. Sanz-Cillero * 

Grup de Fi'sica Teorica and IFAE, Universitat Autonoma de Barcelona, 08193 Barcelona, Spain. 

The use of Pade Approximants for the analysis of the pion vector form-factor is discussed and justified in this 
talk. The method is tested first in a theoretical model and applied then on real experimental data. It is shown 
how the Pade Approximants provide a convenient and reliable framework to incorporate both low and high 
energy information in the euclidean region, leading to improved determinations of the low energy parameters 
such as, e.g., the quadratic radius {r^)v- 



The object of the study of Ref. [T] was the pion 
vector form-factor (VFF), defined by 



(1) 



with the electromagnetic current J'^^^ 
f U7^'u - \dY'd - \lYs and = ~{p' - pf, 
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such that Q > corresponds to space-like data. 

Our knowledge on the VFF from the time-like 
experimental data tells us that the spectral func- 
tion is essentially dominated by the vector reso- 
nance /9(770). The remaining effects are almost 
negligible. On the other hand, it is a well known 
fact that vector meson dominance (VMD) pro- 
vides a good description of the space-like data up 
to relatively high energies. 

The question raised in this work is not so much 
"why this is so" but rather "how can we make 
use of this knowledge" . Our aim is to construct a 
space-like description of the VFF that can be sys- 
tematically improved. The Fade Approximants 
(PAs) provide this suitable framework. 

A Fade Approximant to a given function is 
the ratio of two polynomials Rn{z) and Qm{z) 
(with degree N and M, respectively), constructed 
such that the Taylor expansion around some point 
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(the origin in our case) exactly coincides with 
that of F{z) up to the highest possible order, i.e. 
F{z) - Pi^{z) = ©(z^^+^+i). This may not look 
a big improvement with respect to the usual Tay- 
lor power expansion as they use the same inputs, 
the low energy coefficients. However, the poly- 
nomial is unable to handle singularities such as 
the logarithmic branch cut of the VFF, whereas 
they are partially mimicked by the Faded. Thus, 
in many cases the Fade works far beyond the 
range of applicability of Taylor expansions, al- 
lowing the incorporation of both low and high 
energy data. From this perspective, VMD, where 
F{Q^) = (1 + gV^vMo)"'' is just a Fade po, 
the first term of a sequence of the form Pf". 

Our analysis is focused on euclidean quantities 
in real Quantum Chromodynamics (QCD), i.e., 
with Nc — 3. We separate away here from other 
studies based on QCD in the limit of large number 
of colours |[2j where the large- A^c amplitudes are 
approximated through rational functions [31415) . 
In some situations, the poles and residues were 
real and could be related to hadronic masses and 
couplings. Instead, in the present work the Fades 
are just a useful mathematical tool which does 
not know about particles, masses, widths... Our 
aim is to describe the euclidean momentum range 
and doing that the Fades tend to mimic the shape 



^Notice that if there were no singularities the Taylor ex- 
pansion would converge everywhere in the complex plane 
and the Pades would not provide any further improve- 
ment. 
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of the physical spectral function -strongly peaked 
around the rho mass but finite- through a finite 
number of poles. 

The traditional point of view with respect to 
VMD has been that the rho meson dominates, 
having the remaining resonances small couplings. 
As we said, the PAs do not care about the parti- 
cle poles in the amplitude. The Pade poles rather 
seem to be more related to bumps in the spectral 
function than to hadronic poles (masses, widths, 
couplings...). The danger of relating the poles 
of the rational approximants with mesonic reso- 
nances is clear. 

Obviously, unlike the space-like region, we do 
not expect to reproduce the time-like data since 
a PA contains only isolated poles and cannot re- 
cover a continuous time- like cut. 

A final important point is that the PAs do not 
care about the asymptotic behaviour of the am- 
plitude at very short distances. Although they 
work for a much wider range of momenta than a 
polynomial, the Pade eventually breaks down be- 
yond some energy and they are not intended to 
reproduce the VFF at — > cxo. 

The importance of these techniques is that they 
are able to produce alternative outcomes compet- 
itive with other approaches. The Pades are con- 
structed in a simple and systematic way. They 
provide an efhcient procedure which allows in- 
corporating high-energy space-like information. 
This makes it an interesting alternative tool for 
the analysis of euclidean amplitudes. 

There are several types of sequences of PAs 
that may be considered. In order to achieve a 
fast numerical convergence, the choice of which 
one to use is largely determined by the structure 
of the function to be approximated. In this re- 
gard, a glance at the time-like data of the pion 
form factor makes it obvious that the form fac- 
tor is clearly dominated by the rho meson con- 
tribution. The effect of higher resonance states, 
although present, is much more suppressed. In 
these circumstances the natural choice is a single- 
pole Pade sequence Pf" [6], i.e. the ratio of a 
polynomial of degree L over a polynomial of de- 
gree one. Nonetheless, one should not confuse the 
Pade pole with the rho mass. 

In order to test the aforementioned single- 



pole dominance, we have also considered the se- 
quence P2 , confirming the results found with 
the PAs Pi. We have also considered the so- 
called Pade- Type approximants (PTs) I4I7J (ra- 
tional approximants whose poles are fixed before- 
hand; here they have been taken to be the physi- 
cal masses). Finally, we also considered an inter- 
mediate case, the so-called Partial-Pade approx- 
imants (PPs) 7J, in which some of the poles are 
fixed to the physical masses and some are left free 
and fitted. The different results were rather inde- 
pendent of the kind of rational approximant used, 
being all of them consistent among themselves. 

Test with a model 

In order to illustrate the usefulness of the PAs, 
we will first use a phenomenological model as a 
theoretical laboratory to check our method. 

A VFF phase-shift with the right threshold be- 
haviour is considered ,8.9il0j : 



S{t) = arctan 



M2 



t 



(2) 



with the t-dependent width given by |8l9j 



rp(t) = ro 



^3(M2) 



e{t-^ml) , (3) 



and (j{t) = \/l — Ain^/t. The input parame- 
ters are chosen to be close to their physical val- 
ues Fo = 0.15 GeV, = 0.6 GeV^, Ami = 
0.1 GeV^. The form- factor is then recovered 
through a once-subtracted Omnes relation, 

p 



F{Q') = exp ■ 



w^S(t + Q2) • (4) 



At low energies, the form-factor is given by the 
Taylor expansion. 



1 + ^ak{-Q 



2\k 



(5) 



where the coefficients ak are known since they are 
determined by Fq, Mp and 4toJ. The condition 
F(0) = 1 has been already incorporated. 

In order to recreate the experimental situ- 
ation [TT]-[Tn], an emulation of the space- like 
experimental data is generated in the range 
0.01 GcV2 < Q2 < 10 GeV2 |T]. 
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Figure 1. The sequence of Pi PAs is compared 
to the available space- like data [11]- US]: (brown 
dashed), pI ( green thick-dashed), Pi (orange dot- 
dashed). Pi (blue long-dashed), Pi (red solid). 



The Fade Approximants Pi{Q^) are then fit- 
ted to these euclidean "data" points, providing a 
prediction for the low-energy coefficients a^.. It 
is found that as L increases the sequence of FAs 
Pi converges to the exactly known results, al- 
though in a hierarchical way, i.e. much faster for 
ai than for 02, and this one much faster than 
as, etc... The relative error achieved in deter- 
mining the coefficients ak by the Fade Pi was, 
respectively, 1.5% and 10% for ai and 02. These 
results will be taken as a rough estimate of the 
systematic uncertainties when fitting the real ex- 
perimental data with Fades in next section, and 
they will be added to the final error [1]. 

Experimental data analysis 

All the available experimental data in the 
space- like region have been employed [llj-|16j. 
ranging in momentum from = 0.015 up to 
10 GeV^. 

As discussed in the introduction, the prominent 
role of the rho meson contribution motivates that 
we start with the Pf" Fade sequence. 

The fit to Pi yields a determination for the 
pion VFF and the coefficients a^. Nonetheless, 
according to Ref. [T7] , the VFF is supposed to fall 
like (up to logarithms) for — > 00. This 

means that, for any given value of L, one may 
expect a good fit only up to a finite value of Q^, 
but not for asymptotically large momentum. This 
is clearly seen in Fig. [U where the Fade sequence 
Pi is compared to the space-like data. The Fades 
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Figure 2. ai and 02 Taylor coefficients for the Pi 
PA sequence. 



converge to the real VFF at low and mid energies 
but eventually diverge. 

Fig. [21 shows the evolution of the fit results for 
the Taylor coefficients ai and 02 ■ As one can see, 
after a few Fades they become stable [1]. Thus, 
our best fit is provided by Pj* [1], yielding 

ai = 1.92 ±0.03 GeV^^ , . 

a2 = 3.49 ± 0.26 GeV^"* , 

with a xVdof = 117/90. 

Alternative rational approximations were also 
constructed [1]: two-pole FAs, Fade- Types and 
Fartial-Fades [4]. They all provided compatible 
results. 

Conclusions 

Combining all these outcomes, we obtained the 
results shown in Table [1] For comparison with 
previous analyses, the value of the quadratic ra- 
dius is provided (given by {r'^)y — Gai, with 
our best determination ai — 1.907 ± 0.010^^^ ± 
0.03 GeV-2). 

sys / 

In summary, in this work we have used ratio- 
nal approximations as a tool for fitting the pion 
vector form factor in the euclidean range. Since 
these approximants are capable of going beyond 
the low-energy region, they are rather suitable for 
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Table 1 

Our results for the quadratic radius {r'^)v and second derivative a2 are compared to other determinations [101181 
I19I20I21I22] . Our first error is statistical. The second one is systematic, based on the analysis of the VFF model 
of the previous section [J . 

{r')l (fm") aa (GeV^) 

This work [5 0.445 ± 0.002,tat ± O.OOTsyst 3.30 ± O.OSstat ± 0.33syst 
CGL [HITg] 0.435 ±0.005 

TY [101 0.432 ±0.001 3.84 ±0.02 

BCT [H] 0.437 ±0.016 3.85 ± 0.60 

PP [10] 0.430 ± 0.012 3.79 ±0.04 
Lattice [221 0.418 ± 0.031 



the description of space-like data. They allow to 
extract the low-energy coefficients, improving on 
their determination by incorporating also the eu- 
clidean high-energy data information. As one can 
see in Tabled! the achieved degree of uncertainty 
is shown to be competitive with previous analyses 
existing in the literature |10|18|19|20|21|22] . 

REFERENCES 



1. P. Masjuan et al. [arXiv:0807.4893| [hep-ph]]. 

2. G. 't Hooft, Nucl. Phys. B 72 (1974) 461; 
75 (1974) 461; E. Witten, Nucl. Phys. B 160 
(1979) 57. 

3. S. Peris, Phys. Rev. D 74 (2006) 054013; 
P. Masjuan and S. Peris, JHEP 0705 (2007) 
040. 

4. P. Masjuan and S. Peris, Phys. Lett. B 663 
(2008) 61^ 

5. A. Pich, |arXiv:hep-ph/0205030| . 

6. G.A. Baker and P. Graves-Morris, Pade 
Approximants, Encyclopedia of Mathematics 
and its Applications^ Cambridge Univ. Press. 
1996; chapter 3, sections 3.1 and 3.2; C. 
Bender and S. Orszag, Advanced Mathemat- 
ical Methods for Scientists and Engineers I: 
asymptotic methods and perturbation theory, 
Springer 1999, section 8.6. 

7. C. Brezinski and J. Van Inseghcm, Fade Ap- 
proximations, Handbook of Numerical Anal- 
ysis, P.G. Ciarlet and J.L. Lions (editors). 
North Holland, vol. III. See also, e.g., C. Diaz- 
Mendoza et al, Appl. Num. Math. 53 (2005) 
39 and references therein. 

8. F. Guerrero and A. Pich, Phys. Lett. B 412 
(1997) 382-388. 



9. D. Gomez Dumm et al., Phys. Rev. D 62 

(2000) 054014; J. J. Sanz-Cillero and A. Pich, 
Eur. Phys. J. C 27 (2003) 587. 

10. A. Pich and J. Portoles, Phys. Rev. D 63 

(2001) 093005. 

11. S.R. Amendoha et al. (NA7 Collaboration), 
Nucl. Phys. B 277 (1986) 168. 

12. V. Tadevosyan et al. (JLab F(pi) Collabora- 
tion), Phys. Rev. C 75 (2007) 055205; 

13. T. Horn et al. (JLab F(pi)-2 Collaboration), 
Phys. Rev. Lett. 97 (2006) 192001; T. Horn 
et al. (JLab) arXiv:0707.1794 [nucl-ex]]. 

14. C. N. Brown al, Phys. Rev. D 8 (1973) 
92; C. J. Bebek et al, Phys. Rev. D 9 (1974) 
1229; C. J. Bebek et al, Phys. Rev. D 13 
(1976) 25; C. J. Bebek et al, Phys. Rev. D 
17 (1978) 1693. 

15. P. Brauel et al, Z. Phys. C3, 101 (1979). 

16. Dally et al, Phys. Rev. Lett. 39 (1977) 1176. 

17. G. P. Lepage and S. J. Brodsky, Phys. Lett. B 
87 (1979) 359; Phys. Rev. D 22 (1980) 2157; 
Phys. Rev. D 24 (1981) 1808. 

18. H. Leutwyler, |arXiv:hep-ph/0212324|; 
G. Colangelo, Nucl. Phys. Proc. Suppl. 131 
(2004) 185-191. 

19. G. Colangelo, J. Gasser and H. Leutwyler, 
Nucl. Phys. B 603 (2001) 125. 

20. J.F. de Troconiz and F.J. Yndurain, 

Phys. Rev. D 65 (2002) 093001; Phys. Rev. 
D 71 (2005) 073008. 

21. J. Bijnens, G. Colangelo and P. Talavera, 
JHEP 05 (1998) 014. 

22. P. A. Boyle et al. ' arXiv:0804.3971l [hep-lat]. 



